An upper bound for asymmetrical spinless Salpeter equations 
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A generic upper bound is obtained for the spinless Salpeter equation with two different masses. 
Analytical results are presented for systems relevant for hadronic physics: Coulomb and linear 
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I. INTRODUCTION 

The spinless Salpeter equation (SSE), whose general form is given by (H = c = 1) 



(VP 2 + m? + VP 2 + ^ + ^(r))|^>=M|0) ) (1) 



where M is the mass of the system, is the simplest relativistic eigenvalue equation. It is sometimes denoted semirel- 
ativistic since it is not a covariant formulation. This equation can be considered as a Schrodinger equation with its 
nonrclativistic kinetic part replaced by a relativistic counterpart. More rigorously, it is obtained from the covariant 
Bethe-Salpeter equation [lj with the following approximations: Elimination of any dependences on timelike variables 
and neglect of particle spin degrees of freedom as well as negative energy solutions [2| . It has also been shown that 
a Nambu-Goto string meson with a short-range Coulomb interaction is described by a SSE with a time component 
i ' vector funnel potential, for fixed angular momentum and large radial excitation |3[ . 

Numerous techniques have been developed to solve numerically this equation [4-8]. Though a quite high precision 
can be reached, it is always interesting to obtain bounds on eigenvalues. First, they can be used as verifications 
of the computation. Second, they are generally a simple and fast way to obtain information about the spectra. 
Sometimes, bounds yield analytical results which can give precious details about the solutions of the SSE. With some 
rare exceptions [9|, [lOj], most of the analytical results have been obtained for the symmetrical version of the SSE 

in: Effl 

O" (Wi? 2 + m 2 + V(r)) \<j>) = M\ct>}. (2) 

The one-body (two-body) case is treated with a = 1 (2) . Note that an arbitrary positive value of a can be considered 
to study duality relations between different many-body systems 12111 . 

In this work, we will use the auxiliary field method (AFM) [lOl l2~iTj23j to obtain a generic upper bound on the 
eigenvalues of the general SSE. The computation are presented in section UU after a brief description of the AFM. 
In section IIII1 some analytical results obtained with this bound are given and compared with numerical solutions. 
Concluding remarks are given in section IIVI 



II. UPPER BOUND 



At the origin, the AFM was introduced to get rid of the square root kinetic operator in calculations for semirel- 
ativistic eigenvalue equations J0, . More recently, it has been generalized to treat a greater variety of problems 
[Tol [2l| - |23j . As it is shown in |26j. the AFM has strong connections with the envelope theory H3 • Nevertheless, 

both methods have been introduced from completely different starting points. In particular, the AFM introduces the 
notion of auxiliary fields (also called einbein fields) which is explained below. 
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We will consider here a potential V(r) depending only on the radial distance. The idea is to replace the Hamiltonian 
in ([TJ) by the following one 

H = *±?* + ^! + | + P P(r) + V(I( P )) - P P(I(p)), (3) 
with three auxiliary holds v-y , v 2 and p, and 

" = ^ «a- *-\*), («) 

The kinematics is nonrclativistic, P(r) is the auxiliary potential and the prime denotes the derivative. If we put in 
H the field values which extremize H, 



p 2 + mf and p = K(r), (5) 

then H reduces to the genuine Hamiltonian in ([TJ. The approximation lies in the fact that these auxiliary fields are 
no longer considered as operators but as constants. In this case, an eigenvalue of H is given by 

M{v x ,u 2 ,p) = + v l±^l + V{I{p)) pP{I{p)) + e(p, p), (6) 

where e(p, p) is an eigenvalue of the nonrelativistic Hamiltonian 

H NR =^-+pP(r). (7) 

The AFM approximation for an eigenvalue of the genuine Hamiltonian is given by M(p\, p) for which 

dM(i/i,iy 2 ,p) _ dM(i> 1 ,v 2lP ) 



dvi dp 



0. (8) 



If V(r) oc P(r), the treatment of the potential is trivial and the AFM reduces to the replacement of the semirela- 
tivistic kinetic energy by the best possible nonrelativistic counterpart. An AFM eigenvalue is then an upper bound 
[lpj . In the general case, a function g(x) can be defined by V(x) = g(P(x)). It can then be shown that an AFM 
eigenvalue is an upper bound on the exact energy if g(x) is a concave function (g"{x) < 0). This property has been 
demonstrated in the framework of the envelope theory |F2h14| , but can be applied as well to the AFM [26[ . 

Interesting results can be obtained if the auxiliary potential is a power law, 

P(r) = sgn(» r p with p > -2, (9) 

whose sign is chosen in order that p is always a positive quantity. In this case, a very convenient parameterization of 
P) is given by [H, H3] 

e (M,P)-^(blp) 2/(p+2) ^ • (10) 



2p V M 

Q is a global quantum number which is exactly known in the following cases: 

• p = 2, Q = Q 2 = 2n + l + 3/2; 

• p = 1, Q = Qi = 2(— a„/3) 2 ' 3 , for I — 0, where a n is the (n + l)th zero of the Airy function Ai; 

• p= -1, Q = Q_i =n + l + l. 

In the general case, accurate values of Q can be easily numerically computed. Simple and good analytical approxi- 
mations can be found in (22[. 

After some algebra, constraints ©, with P(r) given by ©, reduce to 

v} - = (b| K(r ) M ) 2/(p+2) Q 2 ^ 2 ', (11) 
r p = (\p\K(r )p)- p /^Q^+ 2 \ (12) 
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where ro is defined by r$ = I(p). From (jlip and (p~2|) . we can deduce that 

Q x 2 



^1 - TO 1 = ^2 - m 2 = ( ^ ) ' ( 13 ) 



It seems then quite natural to define po = Q/r^ in such a way that v\ = \/p\ + m \ an d ^2 = y/Po + m 2- Taking into 
account that -K'(ro) = V' (ro)/(|p|rQ _1 ), (|TT|) and (fT2")l can be written as 

After some algebra, the approximate eigenvalue © with the parameterization (1101) can be written 

M = ui+u 2 + V(ro). (15) 
Finally, the AFM solution is given by the following set of equations 



M= \lpl + m\ + ^p 2 + m 2 + V(r ), (16) 
Po = (17) 

j} 5 + /2 P ° 2 = ro^'(ro). (18) 

The parameter ro can then be interpreted as a mean distance between the particles and po as a mean momentum per 
particle. If g(x), with V(x) = g(sgn(p) x p ), is a concave function or if V(x) oc sgn(p)a; p , then M is an upper bound 
on a true eigenvalue. It is re markable that the only trace of the auxiliary potential is contained in the value of Q. 
With the notation, T(x) = \/ x 2 + m\ + yj x 2 + mf, the system (|16[) -(fT8 |) can be recast under a form similar to the 
one of the system (3.2)-(3.4) in (23[. This shows that (|18[) is the translation into the AFM variables of the generalized 
virial theorem [28|. 




Figure 1: Classical circular motion of the two relativistic particles. 



A semiclassical interpretation of the system (fT6|) - ([T8l) is also possible. Let us assume a classical circular motion for 
the two particles, as illustrated on fig. [TJ In this case, the force Fi acting on particle i is given by 

F i = m il {v i ) d ^ with ^ = ^. (19) 
at at Ti 

Taking into account that both particles are characterized by the same momentum po, This equation becomes 

Fi = , p2 ° -. (20) 
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The rigid rotation constraint, v\jr\ = V2/T2, implies that 



ro = n +r 2 = n^—^- ; , (21) 

P% + m 2 

with i 5^ j. If the force acting on i comes from the potential V(r) generated by j, then Fx = F2 = V'(ro). ([20)l and 
([21]) can be recast onto the form (fl"8|) . and it is obvious than (fT6|) gives the mass of the system. The total orbital 
angular momentum is rip + r 2 po = TqPo. A semiclassical quantification gives thus r Q p = L + 1/2, and we obtain 
a system very similar to (|16[) - (fT5)l . Nevertheless, the AFM produces more general results: Equations are obtained 
not only for a circular motion; The quantum number Q is unambiguously determined by the choice of the power law 
auxiliary potential; A good choice of this potential can allow to obtain an upper bound; An eigenstate of H^r ([7|. 
with Vi = y/pQ + m 2 and p = K (r ) , is an approximation of the genuine eigenstate (29j . 

III. ANALYTICAL RESULTS 

It is easy to obtain numerical upper bounds of the SSE with the system (|16p ~(|18 p . One can wonder if it is possible 
to obtain analytical solutions? Due to the presence of the square roots, it is only conceivable with very special 
conditions. If one mass is vanishing (mi = and rri2 = m) and if the potential is a sum of power-law interactions 
V(r) = Y,i sgn(Aj) a* r A * with a t > 0, fpTT |> - (fl8 ]l reduce to 



Q 2 



Q + ^^ Q =V Mai< ^ (22) 

This transcendental equation can only be easily solved for V(r) = —a/r+br, in such a way that only r§ appears in (|22l) . 
This corresponds to the funnel potential often used in hadronic physics: The sort range part is a strong Coulomb 
interaction stemming from the one-gluon exchange mechanism while the long range part is a linear confinement 
0, S S 0, HH| • In this case, (|22|) is a cubic equation in t-q whose solution is a very complicated expression, unusable 
in practice. So we will only consider the two parts of the funnel potential separately. 

A. Coulomb potential 

We can expect a massive particle linked with a massless particle solely in a heavy-light meson, but then a confining 
potential is always present. So the system treated here with the sole Coulomb potential is quite academic. Nevertheless, 
equations are interesting to treat. With V(r) = —a/r, (|22|) becomes 

Q +-=£!=== a, (23) 
yJQ 2 + m?rl 

which implies that a — Q > 0. In this case, the solution is 



ro = Q^mER. (24) 



The replacement of this value in (fT6")l yields to 



M =^m{ l -m} (25) 

This result was found in [Io| but without the discussion given in the following. As expected from scaling arguments, 
the ratio M/m is only a function of a and of the quantum numbers. This property is also valid for the genuine 
eigenvalue. This last result can be compared with the AFM solution of the symmetrical SSE $2$ with the Coulomb 
potential —a/r, 



M^=aJl-{£). (26) 
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This result can be found in [HI, [HI or can be derived with the AFM 10]. As expected, in both cases, if m is vanishing, 
the only energy scale disappears and the system cannot exist. So, the ultrarelativistic limit (m — > 0) is not relevant. 
But the asymmetrical case is more intriguing. Indeed, the rest energy is not properly defined in this case when the 
interaction disappears, since one of the particle is massless. From equations above, a solution M exists only when the 
following condition is verified: 

\ < Q < a. (27) 

The right part of this inequality corresponds to the cancellation of the binding (M — > m when Q a), and the left 
part corresponds to a collapse (M — > when Q — > a/2). This can be checked on fig. [2] The mean radius ro — > oo when 
Q — > a, and ro — > when Q — > a/2. A sufficiently strong interaction must exist to bind the system, but unphysical 
values of the mass could appear if it is too strong. With the Coulomb potential, the nonrelativistic limit is not defined 
by m — > oo but by a — > (see (|2^|) '). So, (|2T|) implies that this limit is irrelevant for the solution (j2"5"]) . 



mr M 




Figure 2: AFM solution with a global quantum number Q for a particle of mass m interacting with a massless particle via a 
Coulomb potential —a/r. Left: mean radius ro in a/m unit. Right: mass M of the system in m unit. 



The choicep = —1 (P(r) = — 1/r oc V{r)) implies that the AFM solution is an upper bound with Q = Q-i = n+l+1. 
For instance, with a — 1.2, only the ground state (Q-i = 1) exists. The AFM value (|25|) gives M/m = 0.9798, while 
an accurate numerical calculation gives M/m = 0.8454 @,||. 

The upper bound is not very close to the exact value, 
but (|25|) brings information about the (approximate) analytical structure of the solution. 



B. Linear potential 

A massive particle linked with a massless particle by only a linear potential can be considered as a reasonable 
approximation for a heavy-light meson. With a massless particle, the radius of the meson can be quite large, reducing 
the influence of the strong Coulomb potential. With V(r) — br, (|2"2"|) becomes 

Q+ /n2 Q " 2 a =brl (28) 

The solution takes the following form 



Q Q 2 Q 3/2 /4m 2 
y o 2m z 2m i V o 

It can be easily checked that the number under the square root is always positive. The replacement of this value in 
(HHJ) yields to the AFM eigenvalue. We find more convenient to write it as a function of the AFM solution of the 
symmetrical SSE (|2|) with m = and the linear potential b r [ic| 

M (<T) = 2^/abQ. (30) 
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After some algebra, we find 



.1 \/ [m£ ] \ 2 + 32m 2 ±M ( $ 2) , 



B± = \/l\4 2) A ± + 16m 2 , (31) 



M 



V2 (M^ 2) ) 2 A_ + 16 to 2 (2v^ M (2) + £ H 
16 to B— 



It can be checked on fig. [3] that this quite complicated equation is an upper bound of the exact result for well chosen 
auxiliary potentials: P(x) = x and P(x) = x 2 . 



M (GeV) 




m (GeV) 



Figure 3: Masses M with quantum numbers (I = 0, n = 0) and (I = 0, n = 1) for a particle of mass m interacting with a 
massless particle via a linear potential br, with 6 = 0.2 GeV 2 . Solid line: accurate numerical solution 0,01; long-dashed line: 
AFM solution with Q = Qi (P(x) = as); short-dashed line: AFM solution with Q = Q 2 {P(x) = x 2 ). 



For the linear potential, both ultrarelativistic and nonrelativistic limits are relevant. In the ultrarelativistic limit 
(to <C \/&), the dynamics is dominated by the motion of the two massless particles. Indeed, (|31[) reduces to 

MsaM " 2)+ ^S- (32) 

In the nonrelativistic limit (to 3> Vb), the dynamics is dominated by the motion of the massless particle, the 
other one contributing quasi solely with its mass. As a matter of fact, the computation shows that pip reduces 

to (M {2) = y/HMjP) 



D 

8 TO 



M 



M « m + M^ 1} + v o '_ / . (33) 



These two formulas can also be obtained with the AFM perturbation theory developed in (23[. The ratios M/Mq 

given by p"2"j) and (I3"3l coincide when m/M^ sa 0.34, and the relative error of both formulas with M/Mq given by 
JT|) is then only 5.5%. 



IV. CONCLUDING REMARKS 



The spinless Salpeter equation is a semirelativistic equation mainly used to study two-body hadronic systems. 
Though a quite high precision can be reached for its numerical resolution, it is always interesting to have bounds on 
eigenvalues. A rapid estimation of the spectra is then possible, and sometimes analytical results can be obtained. As 
most of the works have been devoted to the study of the spinless Salpeter equation with two equal masses, we propose 
here an upper bound for systems with two different masses. 

We used the auxiliary field method whose concept, strongly connected with the envelope theory, is to replace a 
Hamiltonian H for which analytical solutions are not known by another one H which is solvable and which includes one 
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or more auxiliary real parameters. Provided that the structure of the Hamiltonian H is a nonrelativistic kinematics 
plus a power-law auxiliary potential, an eigenvalue computed by the auxiliary field method is simply the kinetic 
operator evaluated at the mean momentum per particle po plus the potential energy computed at the mean radius rg ■ 
The product r po is equal to a global quantum number Q depending on the auxiliary potential, and the value of r is 
the solution of a transcendental equation. With a good choice of the auxiliary potential, the value of Q is analytically 
known and the eigenvalue obtained is an upper bound of the genuine energy. 

Numerical computations can be easily performed with this bound, but analytical closed formulas can also be 
obtained for some specific systems. Two such solutions are presented in the case where one particle is massless. 
The two potentials studied present some interest for the hadronic physics: The Coulomb potential and the linear 
confinement. Deep insights about the structure of the solutions can then be studied with the upper bounds. 
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